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LONG  TERM  GOALS 

Long  term  goals  are  to  observe  and  model  wave  and  current  boundary  layer  processes  to  determine 
wave  dissipation  and  wave-bed  interactions  in  coastal  and  nearshore  regions  using  novel 
instrumentation  techniques. 

OBJECTIVES 

The  primary  scientific  objective  of  this  project  is  to  measure  the  bottom  dissipation  of  surface  gravity 
waves  as  they  shoal  across  the  continental  shelf  as  a  component  of  the  Shoaling  Waves  Experiment, 
SHOWEX.  Detailed  observations  of  the  bottom  boundary  layer  (BBL),  resolving  the  thin  oscillatory 
wave  boundary  layer  over  a  rippled  sandy  bed,  have  been  made  at  two  sites  with  differing  wave 
forcing,  mean  currents  and  sediment  bed  types,  to  develop  a  spectral  wave  dissipation  model  for  the 
continental  shelf.  At  each  site,  continuous  maps  of  the  small-scale  morphology  have  been  made  in  an 
area  surrounding  the  detailed  bottom  boundary  layer  measurements  to  study  changes  in  the  bed  in 
response  to  wave  forcing,  and  to  relate  the  effects  of  these  morphology  changes  on  turbulence  in  the 
BBL.  A  spectral  dissipation  model  is  being  developed  to  include  effects  of  a  wave-forced  mobile  bed, 
and  parameterizations  for  low  frequency  currents  including  strong  baroclinic  tidal  and  wind-driven 
currents. 

APPROACH 

During  1999,  two  long  term  inner  shelf  observation  sites  were  established  to  measure  wave  and  current 
forcing  and  detailed  measurements  of  BBL  -  bed  responses.  The  Monterey  Inner  Shelf  Observatory 
(MISO,  Stanton  1999)  has  a  shore-cabled  instrument  frame  at  12m  depth  in  the  southern  end  of 
Monterey  Bay  at  a  site  with  ,15mm  mean  diameter  sandy  bed  and  moderate,  primarily  narrow  band, 
long  period  swell  forcing  conditions.  A  second  site  was  established  in  1  lm  of  water  offshore  from  the 
Duck  pier  between  28  September  and  10  December  1999  during  the  main  SHOWEX  observation 
program.  The  Duck  site  has  significantly  finer  sediment  and  forcing  is  generally  shorter  period,  with 
frequent  energetic  wave  forcing,  and  strong,  low  frequency,  alongshore  currents  in  comparison  with 
the  Monterey  site.  A  year-long  wave  forcing  and  bed  map  timeseries  at  the  Monterey  Bay  site  has  been 
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completed,  providing  a  unique  data  set  to  test  ripple  formation  models  over  a  wide  range  of  seasonal 
wave  forcing  conditions. 

At  each  field  site,  dissipation  rates  in  the  bottom  boundary  layer  have  been  estimated  from  O(lcm) 
resolution,  three  component  velocity  profiles  measured  by  a  Bistatic  Coherent  Doppler  Velocity 
Profiler  (BCD VP)  instrument  (Stanton  1996,  2001a,b).  These  instruments  measure  vertical  profiles  of 
velocity  and  sediment  concentration  over  a  60  cm  range  above  the  bed  at  a  20Hz  rate.  These  small- 
scale  measurements  of  the  bottom  boundary  layer  are  extended  through  the  water  column  with  high 
speed  Broad  Band  Acoustic  Doppler  Profilers.  Wave  dissipation  rates  in  the  mean  current  and  wave- 
forced  bottom  boundary  layer  are  being  estimated  by  decomposing  mean,  wave,  and  turbulent 
components  of  the  three  component  velocity  vector  profiles  to  estimate  the  wave  work  term  that 
dissipates  wave  energy  in  the  WBBL.  The  co-located  measurements  of  the  velocity  vector  profiles  and 
sediment  concentration  also  allow  the  sediment  buoyancy  terms  in  the  TKE  balance  in  the  bottom 
boundary  to  be  estimated  when  sediment  suspension  is  occurring. 

As  the  local  sediment  morphology  can  greatly  influence  the  characteristics  of  the  bottom  boundary 
layer,  a  two  axis  scanning  sonar  altimeter  was  developed  and  deployed  to  quantitatively  measure 
finescale  morphology  over  a  2  by  2  m  area  spanning  an  area  below  the  BCDV  profile  measurements, 
with  2cm  horizontal  and  0.25  cm  vertical  resolutions  and  10-15  minute  sampling.  These  local 
morphology  maps  were  extended  at  the  Duck  site  by  qualitative  2D  side-scan  morphology  images  out 
to  a  10m  radius  around  the  instrument  frame  to  identify  larger  scale  morphological  features.  In 
addition,  cross-shelf  side-scan  swaths  of  larger  scale  bed  features  were  measured  with  a  mulitbeam 
survey  sonar  from  the  R/V  Creed  in  collaboration  with  Fred  Dobson  to  estimate  mobile  bed  changes  to 
strong  forcing  events. 

Three  modeling  efforts  are  in  progress  to  test  morphology  evolution  and  wave  dissipation 
parameterizations.  A  full  DNS  model  developed  by  Don  Slinn  simulates  turbulent  flow  over  a 
sinusoidal  bumpy  bed.  The  model  is  being  run  at  NPS  with  laboratory  scale  waves,  and  the  scaled 
velocity  fields  have  been  resampled  in  time  and  space  with  a  modeled  BCDVP  response  to  assess  the 
robustness  of  the  stress,  shear  and  shear  production  method  used  estimate  wave  dissipation  rates  in  the 
oscillatory  boundary  layer.  The  DUNE2D  boundary  layer  model  for  small-scale  morphology, 
(Andersen,  1999),  has  been  used  by  a  graduate  student,  LCDR  Steve  Martin,  to  model  bedform 
evolution  at  the  Duck  site.  The  model  consists  of  inter-linked  modules  for  flow,  sediment  transport  and 
morphology.  The  flow  module  is  based  on  solving  the  Navier-Stokes  equations  in  the  vertical  plane 
with  k-epsilon  turbulence  closure.  This  2D  model  is  significantly  less  expensive  than  the  DNS  to  run 
on  ocean-scale  2D  domains,  and  is  being  used  to  test  simple  morphology  formation  models  using 
forcing  and  mapped  morphology  from  the  Duck  site.  The  objectives  in  our  modeling  efforts  are  to  test 
and  improve  on  existing  ripple  formation  models  based  on  spectral  wave  forcing  and  prescribed  mean 
currents,  and  to  compare  predicted  wave  stress  and  wave  work  rates  with  the  SHOWEX  observation 
timeseries.  The  BBL  dissipation  formulations  will  be  tested  against  in  situ  cross  shelf  directional  wave 
buoy  observations  being  made  by  Herbers  and  O’Rielly  at  the  SHOWEX  site,  and  incorporated  into  a 
shelf  wave  model  being  developed  by  Phd  student  Fabrice  Arduin  and  Tom  Herbers  (Arduin  et  al, 
2000). 
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WORK  COMPLETED 


Analysis  of  data  from  prototype  versions  of  the  BCDV  and  scanning  altimeter  used  during  SandyDuck 
and  the  SHOWEX  observations  have  allowed  techniques  to  be  developed  for  estimating  Reynolds 
stresses  and  dissipation  rate  profiles  across  the  thin  oscillatory  boundary  layer  above  a  sloping  and 
rippled  bed.  Part  of  this  work  has  been  accomplished  with  theoretical  contributions  from  NICOP 
collaborators  Paolo  Blondeux  and  Geovana  Vittori.  (Blondeaux  et  al  1999  and  2000). 

Results  from  Direct  Numerical  Simulations  of  a  randomly  forced  oscillatory  bottom  boundary  layer 
developed  by  SHOWEX  collaborator  Donn  Slinn  are  being  used  to  explore  how  the  finite  sample 
volume,  vertical  profiles  of  velocity  made  by  the  BCDV  can  be  optimally  used  to  measure  the  wave 
dissipation  rates.  High  resolution  model  output  from  both  a  flat  bottom  run  and  a  sinusoidal  rippled 
bed  have  been  “sampled”  using  the  acoustic  response  of  the  BCDV  to  investigate  what  part  of  the 
stress  and  strain  and  shear  field  is  resolved  by  the  instrument,  and  importantly,  what  turbulence  terms 
contributing  to  dissipation  are  not  resolved.  These  detailed  simulations  suggest  that  the  limited 
resolution  oscillatory  boundary  layer  measurements  made  with  the  BCDV  adequately  resolve  the 
dominant  wave-induced  shear  production  terms  across  the  oscillatory  boundary  layer  for  bed 
morphologies  observed  both  at  Duck  and  Monterey  Bay 

A  processing  technique  has  been  developed  to  analyze  the  long  time  series  of  bed  morphology  maps 
made  by  the  X/Y  altimeter.  The  processing  extracts  dominant  ripple  direction,  crest  lengths,  ripple 
height  and  wavelength  statistics  from  each  map.  A  MS  thesis  project  is  in  progress  to  perform  a 
statistical  anaylsis  of  the  bed  response  to  selected  wave  and  current  forcing  parameters  for  the 
Monterey  Bay  MISO  timeseries.  A  similar  analysis  has  been  completed  fro  the  Duck  SHOWEX  site, 
and  a  manuscript  describing  the  bed  response  and  comparison  with  exiting  laboratory-based  models  is 
in  progress. 

RESULTS 

Analysis  of  the  0.6cm  resolution  vertical  profiles  of  three  component  velocity  measured  by  the  BCDV 
have  allowed  the  near-bed  vertical  Reynold’s  stress  profiles  to  be  estimated  while  minimizing  the 
contributions  of  the  strong  wave-induced  oscillatory  flow  to  the  cross-shore  <v’w’>  correlations. 
Under  moderate  to  strong  wave  conditions,  the  stress  profile  is  resolved  across  the  4  to  10cm  thick 
oscillatory  boundary  layer  allowing  a  Tow  wave  bed  stress  to  be  estimated.  The  resulting  wave  work 
rate  compares  favorably  to  the  integrated  wave  shear  production  term  estimated  across  the  oscillatory 
boundary  layer.  Stanton  2001a  describes  the  BCDV  instrument  and  it’s  application  to  measuring 
turbulent  stresses  and  dissipation  in  the  wave  bottom  boundary  layer. 

An  example  of  streamwise  and  vertical  velocity  component  profiles  above  a  sinusoidal  bed  from  a 
Direct  Numerical  Simulation  (DNS)  is  shown  in  Figure  la.  The  model  code,  which  solves  Navier 
Stokes  equations  without  turbulence  closure  parameterizations  down  to  molecular  viscous  scales,  was 
developed  by  Don  Slinn  and  Tom  Pierro  at  the  University  of  Florida.  It  was  run  at  NPS  with  a 
laboratory  scale  10cm  streamwise,  2.5cm  cross-stream  by  10cm  high  domain,  a  20cm  sinusoidal  20 
cms'1  amplitude  forcing  at  the  top  level,  and  a  1.8cm  amplitde  stream- wise  ripple. 
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a  b 

Figure  la.  An  8  second  timeseries  of  steamwise  velocity  ( upper  panel),  vertical  velocity  ( middle 
panel)  and  turbulent  correlation  u’w’(z,t )  (lower  panel)  for  a  DNS  simulation  of  laboratory-scale 
waves  over  a  sinusoidal  rippled  bed,  sampled  halfwave  up  the  ripple,  lb,  right  panel.  A  60  second 
profile  timeseries  of  cross-shore  velocity  ( upper  panel),  vertical  velocity  ( middle  panel)  and 
turbulent  correlation  u’w’(z,t )  sampled  at  the  location  of  the  blue  profile  line  in  figure  2a.  A 
Reynolds  average  of  p<u’w’>  represents  the  turbulent  stress  at  each  measurment  bin  above  the  bed. 

These  observations  were  made  with  the  BCDVSP  at  12m  depth  during  a  developing  storm  on 
yearday  316  during  the  SHOWEX.  The  mean  orbital  displacement  was  0.9m  dominated  by  groupy  8 

second  swells. 
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Figure2a.  A  mesh  grid  map  of  the  bed  morphology  during  the  observations  shown  in  figure  lb.  The 
Z  axis  is  an  arbitary  displacement  below  the  scanned  altimeter  used  to  map  the  bed.  The  vertical  blue 
line  represents  the  position  of  the  BCDVPS  velocity  profile,  while  the  red  vector  represents  the 
principal  axis  of  the  velocity  outside  the  wave  boundary  layer,  and  the  black  vector  (to  the  same 
scale )  is  the  mean  current  direction.  Sediment  suspension  was  occurring  during  the  stronger  cycles 
of  the  wave  groups.  Figure  2b.  A  semilog  axis  plot  of  the  shear  production  profile  above  the  bed 
during  a  10  minute  interval  spanning  the  timeseries  in  figure  lb.  The  wave  boundary  layer  extends 
from  approximately  5cm  down  to  the  bed,  while  increasing  dissipation  rates  in  the  outer  boundary 
layer  are  associated  with  the  0.35ms'1  mean  current  at  this  time. 


The  laboratory-sized  model  was  adopted  to  keep  the  run-time  below  months  on  an  Origin  2000 
workstation.  The  vertical  velocity  profile  timeseries  of  middle  panel  of  the  figure  la  shows  the  down¬ 
stream  signature  of  strong  vertical  velocities  associated  with  turbulent  eddies  shed  from  the  crest  of  the 
single  cysle  ripple  spanning  the  10  streamwise  period  model  domain.  Taking  the  vertical  velocity  as 
representative  of  turbulent  kinetic  energy,  the  modeled  energetic  eddies  were  restricted  to  about  1.5 
times  the  ripple  height.  Figure  lb  shows  similar  field  observations  (with  differing  scales)  of  velocity 
structure  above  the  bed  morphology  shown  in  figure  2a.  The  field  observations  show  significantly 
more  structure  extending  well  above  the  wave  boundary  layer  due  to  the  presence  of  the  significant 
alongshore  current  and  more  complicated  bed  morphology.  The  DNS  modeling  does  however  provide 
a  valuable  method  of  assessing  our  ability  to  measure  near-bed  dissipation  rates  based  on  a  Reynolds 
stress  profile  and  wave-period  shear  profiles. 

Techniques  developed  to  measure  and  characterize  small  scale  bed  morphology  during  SHOWEX  are 
described  in  Stanton  2001c.  These  analysis  techniques  have  allowed  the  role  of  forcing  factors 
including  mean  currents,  orbital  displacement,  wave  directional  spread,  wave  groupiness  and  swell 
wave  asymmetry  on  the  bed  formations  to  be  statistically  inferred  for  the  SHOWEX  duck  site.  These 
morphology  changes  are  now  being  compared  with  selected  timeseries  of  BBL  turbulence  to  assess 
morphology  contribution  to  a  net  friction  factor  which  can  be  applied  to  a  fine  scale  regional  model. 
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IMPACT  /  APPLICATIONS 


Observations  of  cross-shelf  wave  shoaling  and  energy  loss  under  low  wind  conditions  across  the 
continental  shelf  (for  example  Herbers  et  al,  2000)  suggest  that  bottom  dissipation  is  a  zero’th  order 
term  in  the  cross-shelf  wave  evolution.  Modeling  of  bottom  dissipation  in  coastal  regions  will  directly 
improve  shelf  wave  models,  which  have  wide  ranging  navy  and  civilian  applications. 

RELATED  PROJECTS 

This  research  has  benefited  from  and  contributed  to  the  ONR-sponsored  SandyDuck  program  in  the 
development  and  deployment  of  the  scanning  X/Y  altimeter  and  the  prototype  BCDV,  which  address 
overlapping  issues  in  both  programs. 
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